Bohmian trajectories and quantum phase space distributions 
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Abstract 

We prove that most quasi-distributions can be written in a form similar to that of the de 
Broglie-Bohm distribution, except that ordinary products are replaced by some suitable non- 
commutative star product. In doing so, we show that the Hamilton-Jacobi trajectories and 
the concept of "classical pure state" are common features to all phase space formulations of 
quantum mechanics. Furthermore, these results provide an explicit quantization prescription 
for classical distributions. 
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> : 

txjl Introduction 

£ : 

In addition to the standard canonical and path integral quantizations there is an alternative method 



for quantizing classical theories, namely deformation quantization |l]]-[16]. The virtue of this method 
resides in its formal representation which, in many respects, mimics that of classical statistical 
mechanics. A phase-space quasi- distribution H, which describes the state of the system, plays the 
role of the wave function. Instead of operators, observables are represented by ordinary functions 
in phase space. With hindsight, the results hitherto obtained in this context can be summarized 
in a nutshell: quantum mechanics is obtained from classical statistical mechanics by replacing 
the ordinary abelian product of functions by some appropriate non-commutative associative star 
product. To be more precise, the previous quantization prescription permits a complete specification 
of all the relevant objects of phase space quantum mechanics with (thus far) the notorious exception 
of the quasi-distribution. 

The arbitrariness in the choice of star product is related to ordering ambiguities of operators in 
quantum mechanics. We may choose e.g. to symmetrize operator products completely (Weyl order 
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PI), or else to normal order them (creation operators precede annihilation operators). To this myriad 
of ordering prescriptions there is a corresponding number of star products and quasi-distributions. In 
his work, Cohen has sought to construct a general formalism where all possible quasi-distributions 
(and, consequently, all star products) are related by a sort of gauge transformation which leaves 
physically relevant quantities (mean values, probabilities) unscathed. 

Another alternative formulation of quantum mechanics is the de Broglie-Bohm theory which 
displays a set of remarkable properties that are not shared by the quasi-distribution formulation 
[|17]]-[19 namely that it is based on a positive defined distribution and that it allows for a causal 



interpretation in terms of the so-called "Bohmian trajectories". Both the de Broglie-Bohm and 
the quasi-distribution formulation attempt to provide a "classical-like" description of quantum 
mechanics. The two formulations are obviously equivalent as their predictions are identical to 
those of standard operator quantum mechanics. This notwithstanding, it does not diminish the 
conceptual and practical interest of each individual formulation. Firstly, because in some cases it 
may be advantageous to perform a specific calculation in one formulation rather than another, and 
secondly, because they provide conceptual insights in the interpretation of quantum phenomena. 

In this letter we show that some key features of the de Broglie-Bohm interpretation are shared 
by the quasi-distribution formulation, enforcing the idea that these features play a key role in the 
foundations of the quantum theory (previous work studying the relation between the two formula- 
tions, but exploring different points of view can be found in p(J |21], [2^, pl[). We will prove that: 
1) the dynamics of all quasi-distributions (and in particular of the Wigner function) can be cast 
in terms of the de Broglie-Bohm causal trajectories, 2) the classical pure state is the semiclassical 
limit of all pure state quasi-distributions, 3) all quasi-distributions are ^-deformations of a classical 
pure state. 

From a more technical point of view, these results yield new calculation tools that may find 
interesting applications in a wide range of specific problems: 1) they constitute a more intuitive 
and explicit realization of Cohen's classification, providing a straightforward relation between a 
particular star-product and the corresponding quasi-distribution, 2) they formulate the dynamics 
of all quasi-distributions in terms of a system of two second order partial differential equations, 3) 
they provide an alternative implementation of the Wigner- Weyl map for density matrices, 4) they 
yield a (previously missing) quantization prescription for classical distributions. 



2 Quasi-distributions 

We will consider for simplicity one-dimensional systems. Let ip(x, t) be the wave function represent- 
ing the state of the system at any given time t. The wave function is a solution of the Schrodinger 
equation, with Hamiltonian H = ^— + V(x). Let A(x,p) be any observable, with a complete or- 
thonomal set of eigenvectors {|o >} associated with the eigenvalues The average value of A is 

evaluated according to: < A(x,p) >= J dx ip*{x)A (x, —ifi-^j where we assume that (^(x)) 2 

is properly normalized. Moreover the probability for the measurement of the observable A yielding 

*We will consider for simplicity nondegenerate spectra, although this is not crucial. The spectra can be continuous 
or discrete. 
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the eigenvalue a is given by: 



< ip\a > 



dx tjj* (x)lp a ( x ) 



(1) 



where ip a { x ) — < x \ a > is the eigenvector \a > in the position representation. 

Alternatively, the previous postulates of quantum mechanics can be reformulated in terms of 
phase-space quasi-distributions. According to Cohen's classification 0, the /-quasi-distribution 
Ff(x,p,t) for a pure state is given by: 



F f (x,p,t) 



4tt 2 



h 



h 



dddrj dx' p[a/--Ti,t\il>[af + -7i,t\ /(£, V , t )e^'-^, (2) 



where f(£,T),t) is an arbitrary analytic function of £ and rj. For notational convenience, we will 
henceforth omit the time dependence. In terms of F$ the average value of A is given by the very 
appealing formula: 

< A(x,p) >= J dx J dp F^(x,p)Af(x,p), (3) 



where Af(x,p) is a c-number known as the "/-symbol" associated with the operator A: 

h 



Mw) = ^J d ^J dr > Tr{A(x,p)e^ + ^r\Z, V )}e-^- 



trip 



(4) 



The "/-map" Vf : A(H) — ► A(T*M) attributes to any operator A in the quantum algebra A, acting 
on the Hilbert space 7i, the corresponding /-symbol in the "classical" algebra A over the phase 
space T*M. It can be viewed as a Lie algebra isomorphism, provided A(T*M) is endowed with 
a suitable Lie algebraic structure. Indeed, the /-map renders T*M into a non- commutative Lie 
algebra with an associative *-product and a /-bracket according to: 



' V f (A-B)= A f (x,p)* f B f (x,p) 
V f ([A,B})= [A f (x,p),B f (x,p)] f = A f * f B f -B f * f A f 



(5) 



We could equally define a "dual /-symbol" for the operator A: 

n 



V f (A(x,p)) = A f (x,p) = —JdcJd V Tr{A(x,p)e^ + ^f(^ V )}e^ 



x—vqp 



(6) 



In particular, the quasi-distribution F? is proportional to the dual /-symbol of the density matrix: 

= ^rYf(p) = 2hYf^ > >< 0' ^ e can ODV i ous ly define a dual *-product and a dual 
/-bracket, according to: 



Vj(A-B) 



A,B 



A f( x ,p) */ B f (x,p) 
A f (x,p),B f (x,p) 



(7) 



A f *' f B f - B S *'f A f 



We shall call the sector of observables associated with the star product */ and the /-map Vf, the 
observable 's sector. Likewise, the sector of the quasi-distribution, associated with the dual star 
product *'j and the dual /-map Vf will be designated by dual sector. 
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There is an arbitrariness in the choice of function /(£, rj), or otherwise stated, in the choice of 
quasi-distribution and the corresponding /-symbols [p], §. The physically relevant results such 
as eq.(l) remain unaltered. This constitutes the "gauge" invariance mentioned in the introduction. 

It is worth assembling all the formulae concerning a change from a function 77) to another 
function /2(C) v) ( or > from a quasi-distribution F x (x,p) to another quasi-distribution F 2 (x,p)) |[|: 



A x {x,p) *! = /f 1 



(8) 



9a; ' dp 



dp 



A 2 (x,p) * 2 B 2 (x,p) 



A 1 (x,p) *; B x (x,p) = fi /2" 1 *2 

The most celebrated example corresponds to the self-dual case fw(CiV) = 1j a l so known as the 
Wigner-Weyl formulation of quantum mechanics. The *-product (the Groenewold ^-product ||) 
and the corresponding bracket (the Moyal bracket ||) are given by: 



V w (A-B) = A w * 



V w ( 



A, B 



w B w = A w {x,p)e* J B w {x,p), 
[A W ,B W ) M = 2iA w (x,p) sin (fjf) B w (x,p), 



(9) 



where J is the " Poisson" operator: J = J^J^ — ■§--§z, the derivatives d and d acting on Aw and 
Bw, respectively. The corresponding quasi-distribution is the celebrated Wigner function M. We 
shall discuss some of its properties in section 4. 

Finally, to obtain the counterpart of eq.(l), we first need to introduce the *-delta function, which 



is a non-commutative generalization of the ordinary delta function [15|, [T 

5f(A(x,p)) = ^-jdk ef/^). (10) 
The *-exponential is defined by: e^T = J2™=o h [B(x,p)] n * f , where [B(x,p)} n * f is the n-fold /- 



star product of B(x,p). The *-delta function is a solution of the *-genvalue equation []T3], 
i.e: A f (x,p) */ d( (A f (x,p) - a) = 5( (A f (x,p) - a) * f A f (x,p) = aS( (A f (x,p) - a), where A f = 
Vf (^Aj . The stargenvalue equation can be obtained upon application of the /-map to the eigenvalue 

equation: A\a >< a\ = a\a >< a\. In the Wigner-Weyl representation the *-delta function is a 
^.-deformation of the ordinary delta function |Tj| : 

S w 



(A(x,p)) =6{A{x,p)) - jei{x,p)F (A(x,p)) - ^ e 2 (x,p)5"' (A(x,p)) + 0(h*), (11) 



where Bfap) = gg-(S) 2 ande 2 (x )P ) = (f f -2^ f + (f ) 2 . In particular, 
we get for the position and momentum variables: 5™ (x — x ) = 5 (x — x ) and 5f {p — po) = 
S(p-Po). 
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We are now in a position to write a formula for | < ip\a > | 2 in the phase space formulation^: 

V (A = a) = f dx f dp F f (x,p)5{ {A f {x,p) - a) . (12) 



3 De Broglie-Bohm interpretation 



In the De Broglie-Bohm interpretation of quantum mechanics JT7j-[[il|, the wave function is written 
in the form: 

ip(x, t) = R(x, t) exp (is{x, t) V (13) 

where R(x, t) and S(x, t) are some real functions. Substituting this expression in the Schrodinger 
equation we obtain the dynamics of R and S: 



&p 
dt 


+ 


a_ 

dx 


as 
dt 


+ 


1 

2m 



^1 

4m 



1 d 2 V _ 1 J_ ( &P\ 
■2 \ 



V dx 



dx j 



(14) 



0, 



where V{x) 



x 



R 2 (x), is the probability distribution. The first equation is a statement 



of probability conservation, with associated flux • The second equation is interpreted as a 
Hamilton- Jacobi equation. The solution S(x, t) represents an ensemble of trajectories, known as 
Bohmian trajectories, for particles under the influence of a classical potential V(x) and a quantum 
potential: 



Q{x,t) 



h 2 1 d 2 R 
2m R dx 2 



4m 



1 d 2 V _ 1 J_ ffflA 
V'dx 2 " 2T 2 \dx) 



(15) 



Moreover, the momenta p of the particles are subject to the constraint: p = The ensemble of 
Bohmian trajectories stem from the position x at time t = with momentum p = |^(x,0) and 
probability V(x, 0) = R 2 (x, 0). Consequently, an ensemble of Bohmian trajectories can be described 



statistically by the phase-space distribution function |17 

F B (x,p,t) = R 2 {x,t)5 (p- 
its dynamics being governed by the equation: 



as 

dx 



x, t) 



(16) 



dF B { v 2 1 

-Qj-{x,p,t) = ^ + V(x) + Q(x,t),F B (x,p,t)j , 

where {, } p is the Poisson bracket. 

2 Here A denotes generically the observable A(x,p) irrespective of the particular representation, i.e. irrespective 
of the choice of function f (£,,?]). 
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4 *- deformed causal form of quasi-distributions 



In this section we study some common features of the de Broglie-Bohm and the quasi-distribution 
formulations as well as their relation with classical statistical mechanics. For simplicity we shall 
provisionally focus on the self-dual Wigner-Weyl case (fw(£,,Vit) = !)• 
The Wigner function @] associated with ip(x,t) is given by: 



F 



w 



1 



Z7T J- 



+ 00 



dye~ iyp tjj* x hy, t)ip[x + -fiy, t . 



;i7) 



From the Schrodinger equation and eq.(17), it is possible to obtain the differential equation that 
dictates the dynamics of F w : 

of w i 



ot 



ih 



H 



M 



It 



where H w (x,p) = ^-+V(x) is the Weyl symbol of the quantum Hamiltonian. The Wigner function 
is a real function and admits the marginal distributions: 



/+oo r+oo 
dpF w (x,p) = \4>(x)\ 2 , V(p)= dxF w (x,p) = \<P(p)\ 2 , 
-oo J —oo 



(19) 



where <p(p) is the Fourier transform of ^(^OB From eq.(19) one could be tempted to interpret the 
Wigner function as a true probability distribution in phase space. However, this interpretation is 
immediately spoiled, if one realizes that it can take on negative values. Furthermore, and contrary 
to what happens in the de Broglie-Bohm formulation, the dynamics of the Wigner function does 
not, in general, allow for an interpretation in terms of causal trajectories (except for quadratic 
potentials). 

We will now show that the Wigner distribution (17) and the de Broglie-Bohm distribution (16) 
are in fact an ft-deformation of each other. To do this it will prove useful to expand the Wigner 
function in powers of Ti: 



F w (x,p,t) 



J_ V OC V DO (-1)" (h\ n+m d"^*d™<P rr jyp n+m 

2vr ^n=0 Z^ m =0 n !m! \2 J dx n dx m ■> c V 



cxp 



fi(5-i)]W.# I t) 



(20) 



If we substitute tp = Re^ s in the previous equation, we obtain to order Tl 6 



F w (x,p,t)=R 2 (x,t)5^p-^y^ 



dlV 

dx 



- R 



d 2 R 

dx 2 



5 { p -dx-) + 2i R W l 



dx 



-0(h 4 



(21) 

Notice that odd powers of H do not appear in the previous expansion, because the Wigner function 
is real. We conclude that the Wigner function can be regarded as a ^-deformation of the Bohmian 
distribution. If we want the two to match, we have to impose successive corrections to vanish, i.e. 



ftd 2 R _ (9R\ 
dx 2 I dx J 



and 



a 3 s 

dx 3 



0. There are other conclusions that can be drawn from the expansion 



(21) of the Wigner function: 



3 This result is a consequence of the fact that 5^ (x — xq) = S(x — xq) and 5^(p — po) = S(p — p ). 
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(i) If the corrections to the de Broglie-Bohm distribution in eq.(21) and the quantum potential in 
eq.(15) are negligible (this corresponds to a formal h — > limit), we conclude that: 

Fj(x,p) = R 2 d (x; t)5 L - ^{x- t)) , (22) 

where R c i and S c i obey the classical continuity and Hamilton- Jacobi equations, respectively. Con- 
sequently, this is positive defined and obeys the classical Liouville equation. The expression 
on the right-hand side of eq.(22) is usually regarded as the classical limit of a quantum mechanical 
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(ii) In most cases two distinct quasi-distributions can be expressed as ft-deformations of each other. 
Equation (21) could then be a hint that the de Broglie-Bohm distribution (16) is actually just 
another quasi-distribution. Indeed, we proved in a previous work that it can be implemented in 
quantum phase-space with a particular choice of Cohen's /-function [(2"4|. 

(iii) Perhaps it is possible to re-express the ft-expansion (21) in terms of the star product *w 

Indeed, the conjecture stated in (iii) is a particular case of the following, more general theorem: 
Theorem: Consider any of Cohen's functions f(£,,r]), such that: 

/(£,0) = 1, |^(0,0) = 0. (23) 

Then the following formula holds: 

Ff(x,p,t) = R(x,t) *' f 8i (p-^( x ^ */ (24) 
where R(x,t) and S(x,t) are the solutions of the Bohm equations (14). □. 

The proof of the theorem can be found in the appendix. The previous formula constitutes an alter- 
native (and, hopefully, more intuitive) formulation of Cohen's classification of quasi-distributions, 
which can be stated as follows. Consider an arbitrary pair of real and continuous functions R 
and S satisfying the usual requirements of the causal interpretation [pjj. Given an arbitrary non- 
commutative associative ^'-product, such that A(x) *' B(x) = A(x)B(x), then the corresponding 
pure state quasi-distribution associated with ip = Re^ s is given by equation (24). 

Furthermore, this formula provides (i) a reformulation of the /-dual-map for density matrices 
(i.e. given R and S, we automatically have the pure state /-quasi-distribution (24) which is the /- 
dual-map of the corresponding density matrix.), (ii) a previously missing quantization procedure for 
distributions. The statement that quantization is formally just a substitution of the standard prod- 
uct by a non-local star-product has so far only applied to the observable's sector. Now its validity 
has been extended to the dual sector as well. In the de Broglie-Bohm formulation quantum and clas- 
sical pure states differ only in that the Hamilton- Jacobi functions S and S c i appearing in equations 
(16,22) obey Hamilton- Jacobi equations with or without quantum potential, respectively. Likewise, 
any quasi-distribution is a *-deformed quantization of a classical pure state plus a replacement of 
S c i by S. Notice also that the de Broglie-Bohm quantization prescription is perfectly compatible 
with the deformation procedure. In fact, we can recover (16) from (24) for a particular choice of 
/-function and associated *-product f24|. Finally, in this context we also see that the concept of 



classical pure state is a common feature of the semiclassical limit of all quasi-distributions. 
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Moreover, the " quasi-causal" form of the quasi-distributions (24) provides an alternative formu- 
lation of the dynamics in terms of Bohmian trajectories. To obtain the time evolution we just have 
to solve eqs.(14) and substitute the solutions R(x,t) and S(x,t) in eq.(24). In our opinion this is a 
much more elegant approach. In the usual quasi-distribution formalism, time evolution is obtained 
by solving the *-deformed Liouville equation: = ^ H^F-f , which is typically an infinite 

order partial differential equation and displays a set of solutions that do not allow for an interpreta- 
tion in terms of causal trajectories. Equation (24), on the other hand, casts the time evolution of a 
quasi-distribution as the solution of a system of two second order partial differential equations and 
allows for a straightforward interpretation in terms of "interfering" Bohmian trajectories. Moreover, 
expansions like the one in eq.(21), seem to be related to the topic of Wigner trajectories, which 
correspond to order-by-order quantum corrections to the classical Hamilton trajectories. 

Finally, let us make a comment on the non-commutability of star products. Let us start with the 
de Broglie-Bohm distribution (16). It corresponds de facto to an abelian star product in equation 
(24). Under a gauge transformation we obtain say the Wigner function with the non- commutative 
Groenewold star product (9). What have we benefited from this gauge transformation? We have 
removed the wave function from the observable's sector [17], and the price to pay was to introduce 
a non-commutative star-product in the dual sector. So, apparently, the role of non-commutativity 
is to replace the action of the pilot wave on the observable's sector, i.e. on quantities such as kinetic 
energy or angular momentum. 



5 Example 



To illustrate the previous result (24), we consider the simple example of a free gaussian wave packet. 
At t — the wave function reads IT' 



x i 
4og h 



ip (x) = (2tto-^) 4 exp 
After a lapse of time t, we get from the free Schrodinger equation: 

, , . / 2 \~3 (x — Ut) 2 if 1 

ip(x, t) = [2irs t J exp ^— h ^p \ x - -ut 

where St = ao (l + 2 ^. 2 ) and u = ^. The amplitude and the phase function read: 



R(x, t) = (2-ko 2 ) 4 exp 



(x-ut) 2 



4<T 2 



S(x, t) = -§ arctan (^) + p (x - \ut) + ^,(x - ut) 2 , 



(25) 



(26) 



(27) 



where a = \st\. They solve the de Broglie-Bohm system (14). Substituting the previous expression 
for S(x,t) in eq.(24) with / = fw = 1, we obtain: 



F w (x,p,t) = R(x,t) * w s: 



P-Po 



4ma q<j 2 



(x — ut) 



*w R(x, t). 



(28) 
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Given the fact that |^ is (in this case) a linear function of x, we conclude from (11) that the *-delta 
function reduces to the ordinary delta function. Furthermore, if we execute the first ^-product, we 
get: 

F w (x,p, t) = o h. (f T P^ (n) [P-Po- i£^(* - «*)] } *w R(x, t) = 
= {£ E~ h. (f T £# / dy(iy)"exp [iy (p - p - 5^(x - ut))} } * w J2(x, t) = 



{^JdyR (x- §y,t)exp 



«y p-po 



4mffjff 2 



utfj\ I *vy R(x,t). 



A similar procedure for the second *-product yields: 
F w (x, p, t) = i E^Lo £ (-f )" / <foi2 - f y, t) (^) n P exp [iy (p - p - ^,(x - ut 



± J dyR (x - |y, t) i? (x + |y, t) exp [iy (p - p - (x - 



A brief calculation shows that (cf.(27)):| S (x + |y, — 5 (x — |y, i 
Upon substitution of the previous expression in equation (29), we recover (17). 



(29) 



6 Appendix 

To derive eq.(24), let us first prove the ensuing lemma. Let us start by considering the Mehta 
function |Tl[| , 

F s (x,p) = -^*W(?)eH (30) 

corresponding to fs{£, tj) = e~^^ v . This rule is associated with the standard ordering according to 
which all powers of the operator x precede all powers of the operator p [|IJ. Here (p(p) stands for the 
Fourier transform of tp(x): <f>(p) = ^= J dx i(j(x)e~n xp . We now claim that the following lemma 
holds: 

Lemma: The quasi-distribution F s (x,p) can be written in the form: 



dS 



F b (x,p) = R(x) *' s 5> \p-—( x )j *' s R(x), 
whereQ: 

A(x,p) *' s B(x,p) = A(x,p)e lh dpd* B{x,p). u 
Proof of the lemma: Let us start by expressing if)(x) in the polar form (13): 

F s (x,p) = -^-R(x)e-^ ( dx' Rh>)e* s W-fo*-*\ 



(31) 



(32) 



(33) 



4 Notice that *' s is the dual star product of *as, since fAs{£,,v) = e '^'' = fs 1 (^ 7 l)- It corresponds to the 
anti-standard ordering according to which the p's precede the x's. 



9 



On the other hand, if we expand S(x') in powers of (x' — x) we obtain: 



% -{S{x')-p{x'-x)}= % -\s{x) 



dS 



1 d n S 



p _ _ w j + -^.J. (34) 



Consequently: 



F s (x,p) = ^ K R(x)fdy R( x + y)e-% y ( p -^ {x) ) + ^=^^ {x)yn 

= R{x) 4 Uj dy e -^-Hw)+iEr =2 ^r£#w(%r) 4 ^ 



(35) 



where we expanded R(x + y) in powers of ?/ and used the definition of *' s (eq.(32)). It remains to 
prove that the term in curly brackets is indeed equal to: 



OS 



dx 



2tt 



iy\p- 



(36) 



Let A{x,p) = iy(p- §f(:r)). We then have: e^ P ° e ' {x) > = E£°=o^n, where Q = 1 and 
Q n+ i = A *' s Q n = {A — hy-^j £l n . We conclude that: ej 



Hp HW) _ e A-Hy^ Qo = e iyp ( e B+C fio 



where B = —iyjf-(x) and C = —hy^-. Now notice that: 



B,C 



ihy 2 ^4- Moreover, any mul- 



di 



tiple commutator involving 



B,C 



and any non zero number of B's vanishes, e.g. B, C , B 



B , C , C , B 



0. The Baker-Campbel-Hausdorff formula E5J then reduces to: 



,ra=2 



-1^ 



n! 



C 



c 



e B e c , 



(37) 



where for each n, the operator C appears n — 1 times in the multiple commutator. A straightforward 



calculation yields: 



,C ,C = -ih n - l y n ^{x). Consequently: 



e<w , j y (-n y ) n d n s i 
\h^ 2 n\ dx 



(x) e- iy ^ x) e- hy ^VL . 



(38) 



Finally, taking into account that e ^a^f^o = 1, we get: 



dy e ty ( p fcWl+sIii.^' m Sz" W 



(39) 



which is in perfect agreement with equation (35). □ 

As a consistency check, we can verify, using the previous equation, that: p *' s 8*, (p — ff(^)) = 
§§(#) *s 5f/ fp — || (xj). We are now in a position to prove the theorem. 

Proof of the Theorem: Let us start by proving formula (24) for the time-independent case. First 
of all let us verify what are the implications of the constraints (23). For /(£, 0) = 1, it is easy to check 
that Ff(x,p) admits the marginal distribution: / dp F'{x,p) = \if)(x)\ 2 = R 2 (x). Moreover, this 
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constraint entails that any ^-independent function g(x) remains unaltered under a change of function 
v) : f (*Ji> ^) d( x ) = f o) = 9( x )- Finally, the second constraint (23) means that p 
is also invariant under a change of /(£, 77): / i-j^j p = f (o, i^) p = /(0, 0)p + i|£(0, 0) = p. 

The sole purpose of these constraints is to leave §f (a?) and p unchanged under a "gauge" 

transformation. The "gauge" transformation will thus only act on the star product. In particular, 
fw(£,v)i fs(£,v) an d fB(£,v) HU satisfy both constraints. From (8) we then have: 



' dx ' dp 



F s (x,p) 



But, from (8) and (23) we obtain: 



x , 



X 



*' s R(x) 



*8 R (*)} } 



f (i£,i%) {R(x) * w 5f (p - f (*)) * w R(x)} = R(x) *} 5l ( P - 



dx 



X 



*' f R(x) 



(40) 



(41) 



As a consistency check, let us verify that for / = fw, the previous expression coincides up to order h 3 
with eq.(21). From equation (11), we get for A(x,p) = p— || (x): 0i(x,p) = 0, 62(0:, p) 

Consequently: 5^ (p — §§) = 8 (p — 



dx 



~f 24~dxJ° 



dx 

as 



9 3 s , 

cte 3 



p - §f J + £>(ft 4 ). Substituting in eq.(41), we 
do indeed recover (21). Moreover, notice that for S(x) = a + bx + cx 2 , the *-delta function reduces 
to the ordinary delta function. 

To obtain the time-evolution of F*(x,p) (eq.(24)) we just have to follow all the previous steps, 
except that we now use the explicit time dependent polar form of the wavefunction ifj(x, t) = 
R(x,t)exp lS(x,t) .□ 
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